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<N , Abstract 
> . 

CN , We recapitulate a thermodynamically consistent excluded volume hadron 

' gas model and examine its differences with other "thermal models" used in 

I the literature. Preliminary experimental data for particle number ratios in 
the colhsions of Au+Au at the BNL AGS (11^ GeV/c) and Pb+Pb at the 

. . , CERN SPS (160A GeV/c) are analyzed. For equal values of the hadron hard- 

pH I core parameters the excluded volume model gives essentially the ideal gas 



predictions for the particle number ratios, which is similar to other thermal 



H ' models. We observe, however, the systematic excess of experimental pion 

abundances compared to the ideal gas results. This effect can be explained 

in our model by a smaller pion hard-core volume compared to those of other 

hadrons. The absolute values for particle number and energy densities at the 

chemical freezeout are predicted with a simultaneous fit to all these AGS and 

SPS particle number ratios. 
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I. INTRODUCTION 



Preliminary data for nucleus-nucleus (A+A) collisions with truly heavy beams have re- 
cently become available: Au+Au at llA GeV/c at the BNL AGS and Pb+Pb at 160^4 
GeV/c at the CERN SPS A systematic analysis of these data could yield clues to 
whether a short-lived phase with quark and gluon constituents, the quark-gluon plasma, 
exists during the hot and dense stage of these reactions. 

For the studies of matter properties in A+A collisions, it is of vital importance to de- 
termine whether local thermodynamical equilibrium in the system is reached. Assuming 
such a local thermodynamical equilibrium at the final (freeze-out) stage of the process, one 
can calculate the particle number ratios without detailed knowledge of the very complicated 
dynamical evolution history of the system. We remind that the chemical freezeout which de- 
termines the hadron number ratios does not necessarily coincide with the thermal freezeout 
which defines hadron momentum spectra. 

The aim of the present paper is to analyze the preliminary data for particle number ratios 
of Au+Au (AGS) and Pb+Pb (SPS) collisions within the framework of the thermodynamical 
equilibrium hadron gas model. The name "thermal model" has often been used in the 
literature for this type of calculations. We stress, however, that those "thermal models" used 
in the literature to calculate particle number ratios are, in fact, very different. Therefore, 
our first step is to explain our model and clarify its difference from other versions of "thermal 
model" . 

It seems natural to start with the ideal hadron gas at the freeze-out stage. All known 
particles and resonances should be included in this gas and the resonance decay modes to 
the observed particles must be taken into account as well. However, such an ideal gas model 
becomes inadequate in high-energy A+A collisions. The chemical freeze-out parameters, 
temperature T and baryonic chemical potential fib, obtained from fitting the particle number 
ratios at AGS and SPS energies lead to artificially large values of total particle number 
densities which are much higher than any reasonable (or "intuitively expected" ) ones at the 
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freezeout. The total particle number density at the chemical freezeout within the ideal gas 
approach is ritot ~ 4r;,o for the AGS and ritot ~ Stt-o for the SPS, where rio = 0.16 fm~'^ is 
the normal nuclear density. These numbers also exceed any experimental estimate obtained 
from the particle multiplicities and from the measurement of the volume of the system at the 
freezeout with pion-interferometry method. To suppress undesirable large values of particle 
number densities, the van der Waals (VDW) excluded volume procedure is used. We will 
follow in the present paper the thermodynamically consistent excluded volume model of 
Refs. |^,|^ which will be recapitulated in the next Section. Some other "thermal models" 
[4-8] also include VDW "corrections" , but in some ad hoc and inconsistent ways. It is often 
believed that the specific details of these VDW corrections are of minor importance. This 
is, however, not the case. The VDW repulsion does not essentially alter the particle number 
ratios, but it does always give a very strong suppression effect on the values of particle 
number densities themselves. As we shall see below the total values of particle number 
densities are suppressed by a factor of 8 and 14 for AGS and SPS chemical freeze-out states, 
respectively. Therefore, the influence of the VDW excluded volume procedure on particle 
densities, as well as on all other thermodynamical functions of the hadron gas, is very strong, 
and hence the correct form of the excluded volume model formulation should be employed. 



II. VAN DER WAALS REPULSION IN THE GRAND CANONICAL 

FORMULATION 

For non-relativistic statistical mechanics, the use of a grand canonical ensemble is usually 
just a matter of convenience. However, in the hadron gas considered below it is unavoid- 
able. In a relativistic theory one can not fix the number of particles, the number of pions 
in a hadron gas system, for example. The average number of particles only makes sense 
in an equilibrium system, like in the case of a photon gas. This average number is not 
conserved, but increases with increasing temperature of the system. This takes place at all 
temperatures, even at very low ones when the thermal motion of each individual pion can 
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be treated non-relativistically. The canonical ensemble with a fixed number of pions has 
no physical meaning. This statement is also valid, of course, for other hadrons. Particle 
chemical potentials, in general, regulate not particle numbers but the values of conserved 
charges. Strongly interacting matter has three conserved charges, viz., baryonic number B, 
strangeness S and electric charge Q (strangeness is conserved as we neglect "slow" weak 
interactions). Therefore, the canonical ensemble can only be defined with fixed values of B, 
S and Q, not the fixed numbers of pions, nucleons and other hadrons. 

Let's start with the ideal gas of one particle species with temperature T, chemical poten- 
tial n and volume V. The pressure p*"* is related to the grand canonical partition function 
Z^"^ and in thermodynamical limit y — > oo, it has the following expression: 

with 



f{k) 



/(yt2^ ^2)1/2 _ X 

expf^ ^1 +77 



(2) 



being the momentum distribution function, d is the number of particle internal degrees of 
freedom (degeneracy) and m is the mass. The value of is —1 for bosons, +1 for fermions 
and 77 = gives the classical (Boltzmann) approximation. The ideal gas particle number 
density is given by 

Furthermore, the entropy density is defined as s"^ — (9p"^/9T)^. The energy density can be 
found from another thermodynamical relation e^"^ — Ts'^'^ — p^"^ -\- ji 7^*'^ and is written as 

e'\T,p) = I dke{e + mY'm. (4) 

For a fixed particle number N, the VDW excluded volume procedure means the substi- 
tution of the volume of the system V by V — vN, where v is the parameter corresponding 
to the proper volume of the particle. Note that this VDW procedure, interpreted in statis- 
tical mechanics as the gas of hard-sphere particles with radius r, requires that the volume 



parameter v equal to the "hard-core particle volume", |7rr^, multiplied by a factor of 4 0. 
To introduce the excluded volume (a la van der Waals) in our grand canonical ensemble for- 
mulation, we start with the volume substitution of V hj V — vN in the canonical partition 
function Z for each fixed separately. The grand canonical partition function for the ideal 
gas system 

Z^\T, f^,V) = £ exp (^) Z^'^(T, N, V) (5) 

then becomes 

Z{T,^^,V) = f:exp(^) Z'\T,N,V-vN)e{V-vN) . (6) 

There is a difficulty in the evaluation of the sum over in Eq. (^) because of the A^- 
dependence of the "available volume", V — vN. To overcome it, we perform a Laplace 
transform on Eq. and obtain 

Z{T,fi,x) = / dV exp{-xV) Z{T,fi,V) (7) 
Jo 

POO ^ ^ ^ 

= / dV exp(-xK) Z'\T,fi,V) , 
Jo 

where fl = fi—vTx and V = V—vN. The second equality in Eq. (|^) can be easily understood 
as follows. We substitute Z(T,fi, V) by the infinite sum of Eq. (|^), change the integration 
variable V to in each term of the series and then use Eq. (J^) to sum up the series again. 
From the definition of the pressure function 

p(T.^).rii,„!^i^^^ (8) 

one concludes that the grand canonical partition function of the system, in the thermody- 
namical limit, approaches 

p(T,/i) V 



Z{T, jj,, V)\v-^oo ~ exp 



(9) 



From the first equality in Eq. (|^) one sees that this exponentially increasing part of Z{T, fi, V) 
generates an extreme right singularity in the function Z{T,fi,x) at some point x*. For 



X < p/T the integration over V for Z{T,fi,x) diverges at its upper limit. Therefore, the 
extreme right singularity of Z{T,fj,,x) at x*{T,fi) gives us the system pressure, 

p(T,/i) = Tx*(T,/i) . (10) 

Note that this direct connection of the extreme right x-singularity of Z to the asymptotic 
behavior oo of 2 is a general mathematical property of the Laplace transform. Using 
the above equations we find 

x*{T,fi) = Jim ' ^' — ^ , fii = ^ — V Tx*{T,fi) . (11) 

V^oo V 

With Eqs. (|T|) and (|TD|) we find from Eq. (|ri|) the transcendental equation for the pressure 
p{T, fi) of the gas with VDW repulsion in the grand canonical ensemble: 

p(r,/x) = p^'^(T,/i) , /i = ^ - vpiT^fi) . (12) 

At this point some comments are appropriate. First, we remind that the singularity x* of 
the Z{T,n,x) function has nothing to do with phase transition singularities of the system. 
A singularity at x*{T,fi) exists for the ideal gas as well — one can easily recover the ideal 
gas formulae by putting f = in the above equations. Another result which can puzzle a 
reader is the difference between the familiar form of the VDW repulsion, 

p {V - vN) = N T , (13) 



and our Eq. (|T5). According to Eq. ( pTSD the ideal gas pressure, p*'^ = NT/V, increases due to 
the excluded volume repulsion as expected intuitively while Eqs. (1) and ([T^) clearly shows 
that the ideal gas pressure decreases with the inclusion of the excluded volume repulsion. To 
resolve this "paradox" , one should keep in mind the important difference between standard 
VDW gas treatment with Eq. (0) for fixed particle number N and our system with fixed 
chemical potential fi in the grand canonical ensemble. To clarify this point, let us calculate 
the particle number density: 
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To have a direct comparison with Eq. (0) we now make the Boltzmann approximation, 
i.e., with ?7 = in Eq. (^. We shall see later that at the AGS and SPS chemical freeze- 
out, the Boltzmann approximation leads to a very good agreement (only a few percent 
deviations) for all thermodynamical functions with those calculated with exact Rose and 
Fermi distributions. In Boltzmann approximation, the distribution function f{k), as well as 
all thermodynamical functions of the ideal gas, possesses a momentum-independent factor 
exp(/i/T). One can easily find that 

p(r,/i) = exp (- ^^P^\ p^^(T,f^) , p^'{T,i^) = T n'^(T,/i) . (15) 



After some simple algebra, Eqs. (p^ , p!5D give 



p{T, fi) [l~v n{T, /i)] = n{T, fx) T. (16) 



The form of Eq. (16) amazingly coincides with that of Eq. (13), and therefore, we have no 



contradiction with the standard VDW gas physics. The crucial difference is that we have in 



Eq. ([T6| ) a n(T, /i) function instead of a fixed number N (or a fixed particle density N/V) as 
is the case in Eq. (plsl) . In the grand canonical ensemble with fixed temperature and chemical 
potential, the VDW repulsion leads to a strong suppression of the particle number density. 
It is just this suppression on the number of particles which leads to the decrease of the ideal 
gas pressure after the VDW repulsion begins to come into effect. Eqs. (p!6D and (0) have 
similar form, but their physical meaning are very different. Eq. (|TB|) calculates the pressure 
as the function of A^, V, T while Eq. (|TB]) gives particle number density n(T,^), but only 
after Eq. (|r^) for p(T, /i) function is solved. 

With p(T, jj) as the solution of Eq. (|T^ the particle number density is given by Eq. (^^. 
The entropy and energy densities are 

_ f dp{T,fi) \ s'\T, fi) 

= V^)^ = 1 + .n-(T,/i) ' 



Eqs. ([T^ , p!7| , [T8|) reveal two suppression effects on the particle number, entropy and energy 
densities because of the VDW repulsion, namely 

1) . The modification of the chemical potential /i — /i as given in Eq. ([T2|). In Boltzmann 

approximation it just leads to an additional factor exp( — v p/T). 

2) . A suppression factor [1 + t; n*°'(T, /x)]^^ < 1. 

For a ideal gas system of several particle species i = 1, ...,h, the thermodynamical func- 
tions are additive and equal to the sums of their partial values for different particle species: 

p^''{T,fi,,...,fi,) = j2pTiT,fii) , (19) 

i=l 

and similar expressions for s*"^, n^'^ and e*'^. Note that the index i includes all the information 
about the i-th particle, m^, di, rji, fii, etc. 

The extension of the excluded volume procedure for several particle species is straight- 
forward. The grand canonical partition function of the ideal gas equals to the product of 
Zl'^(T, Hi,V) for each particle species "i". The excluded volume grand canonical partition 
function for several particle species, i = 1, h, with proper volumes Vi, can then be 

written as 

Z{T,^^„...,^^,,V)= ^ ... ^ H^W (^) Zf{T, N,,V) e{V) (20) 

Ari=0 Nh=0 i=l V J / 

with available volume V = V — Y.i=i ViNi . The Laplace transform of Eq. ( PP] ) gives 

frj. N _ 1- lnZ(T,/xi,...,/i,„1/) 
p{T,fii,...,fih) = T \im =p (T,/ii, 

(21) 

i=l 

fli = Hi - Vi p{T,iJi,...,jjh) , i=l,...,h. (22) 

Particle number density for the i-th species can be calculated from Eqs. (21,22) and is found 
to be 
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The total particle number density is the sum of the partial values nj from i = 1 to i = h. 
Note that particle number density rii depends on all proper volume parameters vi,...,Vh- 
The partial pressures Pi = pfiT, jli) are reduced to the ideal gas ones pf'(T, /ij) in the limit 
Vi — > 0. This, however, does not take place for Ui because the particle number density rzj 
with Vi = still feels the presence of other particles with Vj ^ due to the suppression 
factor [1 + Ej=iVj n'f{T,iJi)]-' in Eq. (^. 

Moreover, the total entropy and energy densities of the VDW hadron gas are given as 

s(T,,„...,,,) = ^-j^ - (24) 

i=i 1 + ^j=i'"j nf{T,Hj) 



III. COMPARISON WITH OTHER "THERMAL MODELS" 

We recapitulate in last Section a thermodynamically consistent formulation of the VDW 
repulsion in the grand canonical ensemble. The problem appeared to be not trivial and 
many thermal model formulations of the VDW "corrections" used in the literature do not 
meet the requirement of self-consistency. The essential difference between the model that we 
use and those of Refs. [4-7] is the modification of the hadron chemical potentials according 
to Eq. (0). The VDW "correction" of the ideal gas formulae by the suppression factor 
[1 + f j n*'^(T, was postulated in Refs. [4-7] for all thermodynamical functions 
including the pressure. These "corrected" functions, however, do not satisfy fundamental 
thermodynamical relations. Besides, the consistent formulation allows the dependence on 
the excluded volumes fj's for particle number ratios, while in the formulation of Refs. [4-7], 
particle ratios are always the same as in the ideal gas for any choices of Vi values because 
the suppression factor in the denominator is identical for each i-th particle and therefore 
cancel each other out. 

9 



In Ref. 0, the VDW "correction" was postulated in the form of 

The total pressure p is just the sum of above partial pressures over all particle species i. As 
we have shown in the previous Section, the VDW formulation leads to the transcendental 



equations (^,^) for the pressure function. To compare it with ansatz of Eq. (|26|) , let us 
make again the Boltzmann approximation in Eq. (|2ll) . Using Eqs. (5) and one finds from 
Eq. (|lD that 



p(T, /XI,. ..,/..) - E exp( - ^) pr(T,/xO (27) 
Assuming Vi p/T <^ 1 so that exp (— t>i p/T) ^ {1 + Vi p/T)^^, we use then the expansion 

p(r,/.i,...,/iO = i?''(T,/ii,...,/iO + oiv,) 

and the relation 

p^\T,^i^,...,^i,) = j2p^iT,fi,) - T X:nf(T,/i,) (28) 
i=i i=i 

(the second equality in Eq. (^Sj) follows because of the Boltzmann approximation) to find 
finally 



MT^f^.) ^ . , (29) 



1 + T!;=inf{T,ii,)- 
Eq. (p9|) is still rather different from the prescription of Eq. (p6|) used in Ref. §]. In Eq. (p9[) 
Pi = p**^ when Vi = 0. This is generally true as can be clearly seen in the last equality of 
Eq. (PI). However, this does not take place in Eq. (p6|) if some other Vj ^ 0. We conclude 
therefore that the "excluded volume correction" in Eq. ( |26|) [|| has nothing to do with VDW 
excluded volume procedure even in the limit when all fj's are small, i.e., in the first order 
expansion over Vi. 
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IV. PARTICLE NUMBER RATIOS AT AGS AND SPS ENERGIES 



As mentioned before, particle chemical potentials /ij regulate the values of conserved 
charges. For simplicity we neglect the effects of non-zero electrical chemical potential which 
were considered in Ref. ^0\. Electrical chemical potential is responsible, for example, for 



7r+7r~-asymmetry when colliding ions are heavy and therefore have isotopic asymmetry, i.e., 
their number of neutrons is larger than number of protons. These interesting effects are, 
however, not large numerically at AGS and especially at SPS energies. 
The chemical potential of the i-th particle can be written as 

fJ'i = hj^b + Si Us , (30) 

in terms of baryonic chemical potential fib and strange chemical potential fis, where bi and 
Si are the corresponding baryonic number and strangeness of the i-th particle. The hadronic 
gas state is defined by two independent thermodynamical parameters, T and /if,. The strange 
chemical potential fis{T,fib) is determined from the requirement of zero strangeness 

h 

ns(T, fib, fis) = ^SiUiiT^fii) = 0. (31) 

i=l 

In addition to T and fib, the thermodynamical functions also depend on the excluded volume 
parameters Vi because of Eq. (P^j). 

Baryon and meson resonances and their subsequent decays to observed hadrons are of 
great importance for the measured particle number ratios at AGS and SPS energies. All 
known resonance states with mass up to 2 GeV are included in our calculations. Preliminary 
experimental data for particle number ratios in the collisions of Au+Au at the BNL AGS 
(11 A-GeV/c) and Pb+Pb at the CERN SPS (160 A-GeV/c) are analyzed. We use the 
compilation of the experimental data which were presented by J. Stachel at QM'96 (see 
Ref. [|lT| and references therein). 



The enhancement of strange particle production in A+A collisions has attracted spe- 



cial attention as a possible signal of the quark-gluon plasma formation |12|. However, the 
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"thermal models" of Refs. [4-7] which, in fact, gives the same results as the ideal hadron gas 
model have been successful in reproducing many of the ratios with strange hadrons measured 
at the AGS and SPS without any additional parameter to control strangeness abundance. 
Instead, they have a problem with pion abundances — for the SPS data it has been observed 
[4-6] that a simple ideal gas model is unable, within a single set of freeze-out parameters, 
to reproduce simultaneously the strange particles yields and anti-baryon to baryon ratios 
together with ratios where pions are involved. Specifically, experimental pion to nucleon 
ratio and ratios of pions to other hadrons are larger than the ideal gas predictions. As pion 
multiplicity can be related to the entropy of the system, the ideal hadron gas model appears 
unable to account for the large entropy per baryon of the freeze-out system. The deficiency 
of pions in the ideal hadron gas calculations has become the main problem in the theoretical 
"thermal model" in the interpretation of the SPS particle production data. Several mecha- 
nisms [|^J^,|l^ have been proposed to remedy this problem, but no satisfactory answer has 
yet been found. 

Our model procedure is as follows. For given values of T and fib, we solve the system of 
equations, Eqs. ( PT|p2| , pTD , to find p(T,fii,) and fis(T,fih)- The values of T and fit are then 
determined from the 'best fit' to particle number ratios not involving pions. These best fits 
to the Au+Au (AGS) and Pb+Pb (SPS) preliminary data are shown in Fig. 1. The chemical 
freeze-out parameters are found to be T = 140 MeV, fit = 590 MeV for Au+Au AGS, and 
T ^ 185 MeV, /i^ = 270 MeV for Pb+Pb SPS coUisions. 

All particle number densities are calculated from Eq. ( p3D for all known stable particles 
and resonances with mass up to 2.0 GeV. The total production yield of hadron i is then 
proportional to the sum of its thermal density and all possible resonance decay contributions 
to that hadron i: 

nf' = m +nf^ = rii + J^Uj "UO , (32) 

where (y{j,i) is the probability (branching ratio) for resonance j to strongly decay into 
hadron i. 
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In our calculations we will examine the dependence of hadron ratios on particle volume 
parameters f j's. Results shown in Fig. 1 correspond to the same proper volume parameter 
Vi = V for all hadrons. In this case particle number ratios are almost independent on the 
value of V and are the same as in the ideal gas and as in Refs. [4-7]. It occurs because of 
the relation 

(T,/i,) = exp ^^^^(T,^,) . (33) 

Note that the above equality becomes exact in the Boltzmann approximation. In this 
case a common VDW 'denominator', ^ + Y,'j=iVjn''f'{T, jlj), and a common 'numerator', 
exp(— t; p/T), are canceled and Eq. (^) leads to 

ni{T, /ib) ^ nf^jT, ^g) 

n,(T,/ib) nf{T,^l,) ' ^ ^ 

The value of the parameter v is, however, still crucial for the absolute values of particle 
number densities as well as for all other thermodynamical functions of the hadron gas. At 
the same fixed T and /i^, all thermodynamical functions of the hadron gas are smaller than 
in the ideal hadron gas and strongly decrease with increasing v. 

We use quantum statistics in our calculations, but some of our qualitative arguments 
depend on the validity of the Boltzmann approximation. We have checked for all particle 
number ratios that the Boltzmann approximation in the ideal hadron gas (i.e., Vi = 0) gives 
an accurate estimate to the corresponding quantum statistics values: ~ 1-3 % for both AGS 
and SPS chemical freeze-out parameters. In our consideration with Vi > 0, each chemical 
potential is shifted by —Vip and the Boltzmann approximation always becomes even much 
better. 

In the case when not all of the f j's are equal, the hadron volume parameters Vi do influence 
the particle number ratios through the modification of Eq. (^) in the particle chemical 
potentials, required by thermodynamical self-consistency. The effect is quite evident that 
hadrons which take up less space (i.e., smaller values of Vi), and hence smaller influence on 
the excluded volume, have the advantage. The particle number ratios of those small hadrons 
to larger ones increase in comparison with the ideal gas results. 
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To solve the problem with pion multiplicities in the VDW model we now introduce 
different hard-core radiuses: for pions and r for all other hadrons (r > r^^). Such a 
possibility with r^r = was considered in Ref. We remind that the excluded volume 

parameters are = 4 ■ ^ rf . Using Eqs. (|32| ) and (0), we then find 

^ exp(-t;^p/T)<"(T,/i^ = 0) + E.y. exp(-t;p/T)nf (T, tt) 
n*°* exp(-fp/T)nf (r,/ii) + J2j^ieyipi-'vp/T)nf{T,fij)a{j,i) 

(35) 

_ exp(/i;/T)<'^(T,/i^ = 0) + E,y.nf(T,/i,)a(j,7r) 
~ nf{T,Hi) + J2j^inf{T,fij)a{j,7r) 

where 

1 Gtt 

/i; = - t;^) p = — {r^ -rl) p{T,^b; r,r^) . (36) 

In Eq. (p6D we add r and in the arguments of the pressure function to remind the pressure 
dependence on particle hard-core volumes through Eqs. ( pi]) and (p^). Eq. (|35D shows that 
the ratio of pions to any hadron i is changed in comparison to the ideal gas calculations. 
From the second equality of Eq. ( P^D the increase of the thermal pion density appears to be 
due to an effective pion chemical potential /i*. In the Boltzmann approximation, it leads 
just to the additional factor of exp(/i*/T) in the pion number density of the ideal gas. Note 
that all chemical potentials fii are transformed to jli = — Vip and therefore /Xtt = —Vt^p 
becomes negative because fij^ = 0. If is smaller than Vi = v, in the ratios riT^/ni, fi'^ looks 
like a positive pion chemical potential in the ideal gas formalism. From this explanation 
of the origin of /x*, it is clear that there is no restriction on its possible values by the pion 
mass, in contrast to the ideal Bose gas where always fi < m. 

In Figs. 2 and 3 different pion to hadron ratios are shown for Au+Au AGS and Pb-|-Pb 
SPS collisions, respectively. Preliminary experimental values are designated by dotted lines, 
while our model results are represented by solid curves, as functions of /i* . The value yU* = 
corresponds to the ideal gas results. We have already demonstrated that particle number 
ratios for Vi = (the ideal hadron gas) remain the same as those for Vi = v = constant. For 
all ratios in Figs. 2 and 3 we find that the experimental values systematically exceed the 
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ideal gas results (/i* =0). To fit data, one needs /i* > and from Eq. ( ^6]) it means r > r^r. 
Our T and fib values are already fixed both for AGS and SPS from ratios given in Fig. 1. 
Assuming different from = r for other hadrons we obtain no changes in the VDW 
model values for the ratios shown in Fig. 1. As pions have no influence on those ratios, they 
remain the same as the ideal hadron gas results. 

At fixed T and fib the value of /i* is a complicated function of and r. The ratios 
in Figs. 2 and 3 feel not specific values of (r^,r) but /i* value. From the preliminary 
experimental data we find /i* = 100 MeV for Au+Au AGS collisions and /x* = 180 MeV for 
Pb+Pb SPS ones. We stress that the problem with pion deficiency observed in Refs. [5,6] 
for SPS energies looks similar to what we find in both AGS and SPS data analyzed above. 
This deficiency of pions for the preliminary data in Pb+Pb SPS collisions is not so drastic 



as in S+Pb SPS data, where the pion to nucleon ratio is close to 8 |jT5[ and is approximately 
2 times larger than the ideal gas result. To have an agreement with data in Figs. 2 and 3 we 
need only about 30% larger pion to other hadrons ratios than those in the ideal hadron gas. 
With r^r < r (and therefore //* > 0) the thermal pion number density would increase by a 
factor of exp(yU*/T). It is about 2 for Au+Au AGS and about 2.6 for Pb+Pb SPS collisions 
to have an agreement with data. The effect for the pion thermal density is therefore quite 
strong. However, it does not strongly alter the pion to hadron ratios, Eq. (|35| ) because at 
both AGS and SPS energies the pion production is essentially dominated by the resonance 
decay contributions. Let us also remind the possibility of the chemical non-equilibrium 



effects for pions discussed in Ref. |]T6[. It would lead to the chemical potential /Xtt > with 
its values being always smaller than the pion mass m^r. For r^^ < r, we have obtained in the 
VDW model the effective pion chemical potential /x* whose value is not restricted by the 
pion mass, and no chemical non-equilibrium effects are required. 
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V. PARTICLE HARD-CORE RADII 



Our fits to the preliminary data of Au+Au (AGS) and Pb+Pb (SPS) on tlie particle 
number ratios with VDW hadron gas model are shown in Figs. 1-3. From fitting the data 
we have found the following model parameters: 

AGS : T = 140 MeV , = 590 MeV , /i; ^ 100 MeV , (37) 

SPS : T ^ 185 MeV , /Zb = 270 MeV , /z^ = 180 MeV . (38) 

The set of all possible values of and r which give the same value of fi* in Eq. (|36|) for 
the AGS or SPS defines a curve in (r^^, r)-plane. The two curves corresponding to the AGS 
and SPS respectively are shown in Fig. 4, together with their intersection point (r^ = 0.62 
fm, r = 0.8 fm). This intersection point is a solution for and r in the VDW hadron gas 
model to fit simultaneously the AGS and SPS data for all particle number ratios. 

In Tables I and II, we show the values of the total meson number density n^, baryon 
number density n^, total particle number density ritot (antibaryons included), total energy 
density e and total pion number density n^°^, which includes both thermal pions and con- 
tributions from resonance decays, at chemical freezeout. The ideal hadron gas (r^^ = r = 0) 
results are given in the first row, and those from our calculations with different (r^, r) values 
along the curves in Fig. 4 are shown in the remaining 4 rows. We emphasize that the last 
rows in the Tables give our results for the intersection point (r^ = 0.62 fm, r = 0.8 fm), and 
that these values are our predictions within the VDW hadron gas model for the chemical 
freezeout at Au+Au AGS and Pb+Pb SPS collisions, respectively. We believe that these 
r7r,r values, along with the corresponding densities, give a reasonable physical solution for 
the chemical freeze-out state in A+A collisions considered. 

VI. SUMMARY 

A self-consistent hadron gas model with the VDW excluded volume is considered and 
critically compared with other "thermal models" used in the literature. This approach is 
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then adopted to analyze the prehminary data of Au+Au (AGS) and Pb+Pb (SPS). Within 
the VDW model, the obtained values of particle number and energy densities, and other 
thermodynamical functions at chemical freezeout are very different from those obtained in 
the ideal hadron gas, as seen in the first and last rows in Tables I and II. Because of their 
strong effects on hadron thermodynamical functions, the VDW gas formulation should be 
properly treated. 

The preliminary data of Au+Au (AGS) and Pb+Pb (SPS) for the particle number ratios 
can be fitted in the VDW hadron gas as shown in Figs. 1-3. The model parameters, given 
in Eqs. ( p7|^8D , lead to the enhancement of pions in pion to hadron ratios as compared to 



the ideal hadron gas model predictions. This enhancement in pions is regulated by /i* (> 0) 
and is explained in the VDW model by a smaller pion "hard-core radius" than those of all 
other hadrons. 

The obtained parameters for Au+Au (AGS) and Pb+Pb (SPS) define two curves in the 
(r^, r)-plane, respectively. These two curves shown in Fig. 4 intersect at the point (r^ = 0.62 
fm, r = 0.8 fm). It is the solution in the VDW hadron gas model to fit simultaneously 
the AGS and SPS data for all particle number ratios. The absolute values of the particle 
number and energy densities in the VDW hadron gas for this solution are listed in the last 
rows of Tables I and II. We stress that these values are much smaller than those in the ideal 
hadron gas at the same T and fib shown in the first rows of Tables I and II. There is an 
experimental estimate for the freeze-out pion number density in Si+Pb central collisions at 



the AGS energy: = 0.063 fm^'^ l^. Our result <°* ^ 0.054 fm^^ for Au+Au AGS 



collisions shown in the last row of the Table I is quite close to this experimental estimate. 

The next step is naturally to fit hadron momentum spectra. It requires the inclusion of 
the longitudinal and transverse collective flow effects. The temperature T = 185 MeV looks 
a little too high for use in the fitting of the transverse momentum spectra in Pb+Pb (SPS) 
collisions. For example, the pion inverse slope parameter is near 190 MeV [0] and the freeze- 



out temperature T = 185 MeV seems to leave almost no room for the transverse collective 
motion effects. However, two facts should be taken into account. First, the temperature 
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determined from particle number ratios is for the chemical freezeout which could be higher 
than the thermal freeze-out temperature used in particle spectra calculations. Second, we 
remind again that there are large resonance decay contributions to pion production. Even 
after the enhancement of thermal pions with //* = 180 MeV, resonance decays contribute 
more than 60% to the final pions. These resonance decays are known to lead to a lower pion 
'effective temperature' (the inverse slope parameter) at the transverse pion mass less than 
1 GeV. 
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TABLES 

TABLE I. Meson number density rim, baryon number density n^, total particle number den- 
sity ntot, energy density e and total pion number density n^"*, which includes thermal pions and 
contributions from decays, at the freezeout for AGS Au+Au collision at 11 A-GeV/c. 



(r^r, r) [fm] 


nm [fm-^] 


nb [fm-3] 


ntot [fm~^] 


e [GeV/fm3] 






0.200 


0.402 


0.603 


0.722 


0.400 


(0.00, 0.50) 


0.097 


0.130 


0.227 


0.248 


0.161 


(0.20, 0.52) 


0.089 


0.120 


0.209 


0.228 


0.149 


(0.40, 0.61) 


0.061 


0.084 


0.145 


0.159 


0.103 


(0.62. 0.80) 


0.032 


0.043 


0.075 


0.083 


0.054 


TABLE II. 


Meson number density n^, 


baryon number density nb, total particle number 


density ntot, energy density s 


and total pion number density n^°^, 


which includes thermal pions 


and contributions from decays, at the freezeout for SPS Pb+Pb collision at 160 A-GeV/c. 


(r^r, r) [fm] 


nm [fm-3] 


nb [fm ^] 


ntot [fm-3] 


e [GeV/fm3] 


^tot 


(0.00, 0.00) 


0.771 


0.404 


1.251 


1.585 


1.213 


(0.00, 0.46) 


0.278 


0.095 


0.391 


0.429 


0.382 


(0.20, 0.48) 


0.244 


0.084 


0.344 


0.378 


0.336 


(0.40, 0.59) 


0.139 


0.049 


0.197 


0.219 


0.193 


(0.62, 0.80) 


0.064 


0.022 


0.090 


0.100 


0.088 
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FIGURES 



Fig. 1: Points are the preliminary experimental data for the particle number ratios (see 



Ref. |Tl| and references therein) for Au+Au AGS and Pb+Pb SPS collisions (in the lower 
and upper part of the figure respectively). The short horizontal lines are the model fit with 
T = 140 MeV, jjb = 590 MeV (AGS), and T = 185 MeV, /i^ ^ 270 MeV (SPS). 

Fig. 2: Pion to hadron ratios for Au+Au AGS collisions. The experimental data (see 



Ref. |TT| and references therein) are shown by the dotted horizontal lines. The solid lines 
are the VDW hadron gas model results as the functions of /x* (see the text for details). An 
agreement with data corresponds to /i* = 100 MeV. 

Fig. 3: Pion to hadron ratios for Pb+Pb SPS collisions. The experimental data (see 



Ref. |TT| and references therein) are shown by the dotted horizontal lines. The solid lines 
are the VDW hadron gas model results as the functions of /x* (see the text for details). An 
agreement with data corresponds to yU* = 180 MeV. 

Fig. 4: The solutions of the equations /i* = constant for AGS parameters, Eq. (P7D, 
(the dashed curve) and SPS ones, Eq. (^Sj), (the dotted curve) shown in (r^^, r)-plane. The 
intersection point is approximately at (r^ = 0.62 fm, r = 0.8 fm). 
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